Kinetic theory of surface plasmon polariton in semiconductor nanowires 
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Based on the semiclassical model Hamiltonian of the surface plasmon polariton and the nonequi- 
librium Green-function approach, we present a microscopic kinetic theory to study the influence of 
the electron scattering on the dynamics of the surface plasmon polariton in semiconductor nanowires. 
The damping of the surface plasmon polariton originates from the resonant absorption by the elec- 
trons (Landau damping), and the corresponding damping exhibits size-dependent oscillations and 
distinct temperature dependence without any scattering. The scattering influences the damping by 
introducing a broadening and a shifting to the resonance. To demonstrate this, we investigate the 
damping of the surface plasmon polariton in InAs nanowires in the presence of the electron-impurity, 
electron-phonon and electron-electron Coulomb scatterings. The main effect of the electron-impurity 
and electron-phonon scatterings is to introduce a broadening, whereas the electron-electron Coulomb 
scattering can not only cause a broadening, but also introduce a shifting to the resonance. For InAs 
nanowires under investigation, the broadening due to the electron-phonon scattering dominates. As 
a result, the scattering has a pronounced influence on the damping of the surface plasmon polariton: 
The size-dependent oscillations are smeared out and the temperature dependence is also suppressed 
in the presence of the scattering. These results demonstrate the important role of the scattering on 
the surface plasmon polariton damping in semiconductor nanowires. 

PACS numbers: 73.20.Mf, 73.22.Lp, 72.30.+q, 71.10.-w 
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I. INTRODUCTION 

Since the pioneering theoretical work by Ritchie 1 - and 
the electron-loss spectroscopy measurements by Pow- 
ell and Swan^ physics of surface plasmon polariton 
(SPP) has been extensively studied for more than five 
decades. 3 ~- SPPs are electromagnetic (EM) waves cou- 
pled to the collective excitations of electrons on the sur- 
face of a conductor ! 10 ! 11 In this coupling, the electrons 
oscillate collectively in resonance with the EM waves and 
hence trap the EM waves on the surface. The resonant 
coupling leads to the SPPs and gives rise to their unique 
properties, such as the enhancement of the surface elec- 
tric fields and the slowing down of the group velocity of 
the EM w&ves~2r— Applications exploiting these prop- 
erties have been widely studied in biosensing^i solar 
cells, 22 quantum information processing. 2 '' 25 subwave- 
lcngth optical imaging and waveguiding devices j 11 ' 26 " — 

The SPPs often suffer from dissipative losses^ 2 - Over- 
coming the losses is crucial for the improvement of perfor- 
mance of many SPP-based devices, such as the fidelity of 
the waveguide and the sensitivity of the single-molecule 
sensor. 21 ' 33 An effective modulation of the losses is also 
highly desirable for active plasmonic devices proposed in 



recent years 



34-36 



Thus, a thorough understanding of the 



damping processes responsible for the dissipative losses is 
essential. Since the SPPs in metals have been the major 
focus for many decades, previous studies on the damping 
processes have traditionally been focused on metals. It is 
found that the dominant damping process is the decay of 
the SPPs into electrons, i.e., the Landau dampingj^Zr— 
The interband transitions and the many-body exchange- 
correlations can have pronounced influences on the Lan- 



dau damping. These influences can be incorporated into 
microscopic models based on time-dependent density- 
functional theory or semiclassical model Hamiltonian of 
the SPP4 2 -"— Calculations based on these models have 
shown good agreement with experiments 

In recent years, doped semiconductors, such as SiC, 
GaAs, InAs, CU2S and Cu2Se, are suggested as promis- 
ing candidates to replace metals in SPP applications^"^ 2 - 
The SPPs in doped semiconductors are characterized by 
their substantial low losses and tunable frequencies 
They are also easier to be manipulated via doping or ex- 
ternal electric/magnetic fields and to be integrated into 
complex, functional circuits^ - — The further develop- 
ment of the SPPs in doped semiconductors requires a bet- 
ter understanding of their damping processes. However, 
the physics involved in doped semiconductors can be 
quite different from that in metals. For doped semicon- 
ductors, although the Landau damping is still believed to 
be the leading damping process at large wavevectors , 60 ' 61 
the charge depletion/accumulation laye r 62 ' 63 and the 
electron scatterin g 64 ' 65 are found to have important in- 
fluence on the damping. Of particular importance is the 
effect of the electron scattering, since the typical scatter- 
ing rate can be comparable to the SPP frequency in semi- 
conductors. However, to the best of our knowledge, this 
effect has only been discussed by using phenomenologi- 
cal relaxation times. 64 ' 65 A microscopic theory exploiting 
this effect has not been established yet. 

In this paper, by combining the semiclassical model 
Hamiltonian of the SPP— ~— and the nonequilibrium 
Green-function approach ) 66 ' 67 we present a microscopic 
kinetic theory to study the damping of the SPP in 
doped semiconductors, within which the relevant elec- 



tron scatterings are treated fully microscopically. The 
main purpose of this work is to understand the influ- 
ence of these scatterings on the Landau damping of the 
SPP. To demonstrate this, we focus here on the SPPs in 
InAs nanowire a 23 i 28 ' 29 i 51 and concentrate on the electron- 
impurity (ei) , electron-phonon (ep) and electron-electron 
(ee) Coulomb scatterings. We find that the scattering 
can have pronounced influences on the Landau damp- 
ing of the SPP by modulating the resonance between the 
electrons and the SPPs. Different scattering has different 
effect on the resonance. The main effect of the ei and ep 
scatterings is to introduce a broadening to the resonance, 
whereas the ee scattering can not only case a broaden- 
ing, but also introduce a shifting to the resonance. For 
InAs nanowires, the ep-scattering-induced broadening is 
found to be the dominant effect. These effects can lead 
to a pronounced influence on the damping of the SPP, 
which can be seen from both the size and temperature 
dependence of the SPP damping: (I) The size-dependent 
oscillations of the SPP damping are smeared out, and 
(2) the temperature dependence of the SPP damping is 
suppressed by the scattering. These results demonstrate 
the important role of the electron scattering on the SPP 
damping. 

This paper is organized as follows. In SecHH we intro- 
duce the semiclassical model Hamiltonian of the SPP for 
semiconductor nanowires and briefly outline the deriva- 
tion of the kinetic equations. We also present an analytic 
solution for the SPP damping which provides a simple 
and physically transparent picture to understand the in- 
fluence of the scattering on the SPP damping process. In 
Sec. IIII1 we discuss in detail the influence of the electron 
scattering on the SPP damping by numerical solving the 
kinetic equations. The importance of the scattering is 
demonstrated by studying its influences on the tempera- 
ture and size dependence of the SPP damping. The an- 
alytic solution is also compared with the numerical ones 
in this section. We summarize and discuss in Sec. IIV1 



II. MODEL AND FORMALISM 

A. Semiclassical model for SPP-electron system 

We consider an n-type free-standing cylindrical 
nanowire with radius R as illustrated in Fig. QJa). The 
z-axis is chosen to be along the wire. Following the semi- 
classical approach developed in previous works,— ~— we 
decompose the Hamiltonian into 



H — Hspp + H c \ + i/spp-ei, 



(f) 



where Hgpp, H c \ and -ffspp-ci arc Hamiltonians for the 
SPP, electrons and the SPP-electron coupling, respec- 
tively. Here we only present the Hamiltonian, leaving 
the details to Appendix |A"1 

For nanowires, there exists one fundamental SPP 
mode with axial symmetry, which has no cutoff at 
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FIG. 1: (Color online) (a) Schematic of the structure of the 
electron field E and corresponding surface charge of the axial 
symmetric SPP mode in a cylindrical nanowire. ej° (q) is 
the dielectric constant inside (outside) the nanowire. The 
red curves with arrows represent the electric fields of the SPP 
mode, (b) The dispersive relation of the axial symmetric SPP 
mode for nanowires with different electron density no and 
wire radius R. The thin grey line marks the energy of the LO 
phonon. 



low frequency and has been found to be important for 
both the terahertz emission and quantum subwavelength 



optics 



2:5.28,29.51 



In this paper, we focus on the dynamics 



of this SPP mode. The corresponding Hamiltonian takes 
the form 



#spp = ^2^g b l b g, 



(2) 



where represents the annihilation(creation) boson 

operator for the SPP, with Q q being the corresponding 
dispersive relation illustrated in Fig. [IJb). Note that we 
set H = 1 throughout this paper. 

The Hamiltonian of electrons can be written as 



H e \ — #°i 



ep 



He, 



(3) 



where the free-electron Hamiltonian is modelled by a 
mean-field potential. The interaction Hamiltonians H e i, 
H ep and H ee represent the ei, ep and ee interaction, re- 
spectively. 

By choosing the mean-field potential to be an infinite 
cylindrical potential well with radius R, the free-electron 



Hamiltonian can be written as 



nka 



£ k C nkcr C nka^ 



(4) 



3 



in which e k = 



is the eigenenergy with m* 
representing the electron effective mass. The composed 
index n = (m, n) labels the electron subband with rh and 
n representing the angular and radial quantum numbers, 
respectively. A™ denotes the n-th zero of the Bessel func- 
tion of the first kind Jm{x). The corresponding eigen- 
states read 



ip nk (p,cp,z) 



JmjXjZ Pi R) imtfi ikz 
/^RJrh +1 (Xf) 



(5) 



The ei interaction Hamiltonian can be written as 

^ = EEc'ftW^^- ( 6 ) 

i kqa 

with Ni being the total impurity number and Pi(q) — 
e ~iqzi jj ere we have assumed that the impurities are 
distributed on the surface of the nanowire with an axially 



symmetric distribution. 



is the matrix element for 



the ei interaction. The ep interaction Hamiltonian can 
be written as 



tfc P = £ £ 

Qq nn'ka 



M, 



Qq [ a Qq I" 



J C llk(j C n'k-qtTl7) 



where cLQq(a,Q q ) represents the annihilation(creation) op- 
erator for the LO phonons, with Q and q representing 
the components of the phonon momentum perpendic- 
ular and parallel to the nanowire. Here we use bulk 
phonons in the present investigation, which is valid for 
nanowires with large diameters! 68 ' 69 Mq™ is the matrix 
element for the ep interaction. It is noticed that although 
surface-optical (SO) phonons 69 also exist in nanowires, 
they are of marginal importance since the correspond- 
ing electron-SO-phonon interaction is rather weak com- 
pared with the electron-LO-phonon interaction for the 
nanowires we consider here. The influence of the SO 
phonons will be further addressed in Sec. IIIIBI The ee 
interaction Hamiltonian can be written as 

Hoc = E E £ C nka C li'k'a' C n'k'+q<J' c nk-qa,(8) 
kk' q nn f erer' 

where V™ n is the matrix element for the ee interaction. 
The SPP-electron coupling Hamiltonian -ffspp-oi can be 
written as 

ffspp-d -EE 9k-k' (bk-k> + &£,_*) 4^*., (9) 

nkcr n' k' 

where g™ is the SPP-electron coupling matrix element. 
In these equations, matrix elements v™ 1 , Mq™ , V™ n 
and g" n are given in detail in Appendix 



B. Kinetic equations 

In this section, we briefly outline the derivation of the 
kinetic equations for the SPP-electron system. The de- 
tails can be found in Appendix [B] 



The damping of the SPP is obtained by studying 
the temporal evolution of a coherent SPP wave packet 
with central wavevector Q s , which can be expressed as 
\ B s) =J2 q P? se ~^ Bql2 e b "^\0) (Ref.[nJ- The line-shape 
function of the wave packet p® 3 is chosen to be p® s = 
sin[{Qs-q)L^2] ^ w j iere £ j g ^ wave p ac ket length. Such 

wave packet is typical in a Fabry-Perot SPP resonator, 
which has been observed in various SPP systems! 25 ' 71 " — 
The amplitude of the wave packet can be described by 
B s = YlqPq s Rq- The kinetic equation of B s is obtained 
from the Heisenberg equation of the SPP annihilation 
operator b qi which has the form 



d t B s (t) 



E 

nn r ,kk' ,er 



P%i k ,9^ k ,G<(n'k',nk;tt), (10) 



the "lesser" electron 
G<(nk,n'k';tt') 



where G<(nk,n'k';tt') is 
Green function defined as 

i{clk>vW) c nk*{t)) (Ref.®. 

The kinetic equation of the electron Green function 
G^(nk,n'k';tt) is derived by using the nonequilibrium 
Green-function approach, which can be written as 



id t - {ei - el) G<(nk, n'k'; tt) 

= £(!?_, + B\) [gfG< {nk, nk' - q; tt) 

nq 

-gfG<(nk + q,n'k';U)\ +I^ n , k ,(t), (11) 



where the first term in the right hand side of the equa- 
tions is the coherent driving term of the SPP, while the 
second term is the scattering term consisting the ei, ep 
and ee scatterings. 

Within the rotating wave approximation relative to 
the SPP central frequency f2 S ) 66 ' 67 we obtain the kinetic 
equations for the SPP-electron system 

KB.® =-i E [Pa(nk,n'k';t)]\ (12) 

nn'a kk' 

d t Pa(nk,n'k';t) = iaj 7 ^ P a (nk,n'k';t) 

+ Wk-k'Pk-k'Bs fna(k) — fn'a(k') + Ink,n' k' > (13) 

with the detuning 



J kk' 



p n — O 

e k s l s 



(14) 



In the above equations, B s (t) = B s (t)e lQ ' st and 
P a (nk,n'k';t) = —iG<{nk,n'k';tt)e lVLat represent the 
SPP amplitude and electron polarization, respectively. 
f n(T (k) represents the equilibrium electron distribution 
which is conventionally chosen to be the spin-unpolarized 
Fermi-Dirac distribution at temperature T. 

It should be emphasized that in the derivation, we 
have treated the SPP-electron coupling g™ perturba- 
tively and linearized the equations by keeping only terms 
up to the linear order of <?"" . Thus the SPP couples 



4 



only to the electron polarization corresponding to the 
off-diagonal electron Green function with respect to the 
electron momentum k and subband index n. 

The scattering term within the rotating wave approx- 
imation can be expressed as I = I cl + I cp + I cc , where 
7 el , I cp and I ee are contributions from the ei, ep and ee 
scatterings, respectively. Their expressions can be found 
in Appendix [Bl 



C. Landau damping process 



In the kinetic equations, Eq. (|T3|) describes the res- 
onant excitation of the electron polarizations, while 
Eq. (IT2t describes the back-action of the polarizations to 
the SPP. The summation in Eq. (|12p implies that even 
without any scattering, the phase-mixing between polar- 
izations with different frequencies can also lead to the 
damping of the SPP, which is the origin of the Landau 
damping. 

To further clarify the Landau damping process de- 
scribed by the kinetic equations, we solve the equations 
without the scattering term I. From Eq. (|T3| . the corre- 
sponding polarization can be solved as 

P a (nk,n'k';t) = -5 k ,_ k _ q gf pf B\(t) 

X IfnAV ~ fn'Ak'W^*' +l ° +)t ~ + i0 + ).(15) 

Note that in the derivation, we have assumed that the 
SPP amplitude B s varies slowly compared to the polar- 
ization P. 

By substituting Eq. (f!5)) into Eq. (H~2l and taking the 
long time limit t — > oo, one gets 

-l 



d t B s /B s = -(r 



(16) 



where u> s represents the frequency shifting, expressed as 

\pfgT'?[Mk')-f n ,{k)]- 1 • 



us = 

nk,n' k' ,qa 



J kk' 



-k—qi 



(17) 



in which the summation is understood as a principal 
value integral. r _1 is the damping rate of the SPP, which 
has the form 



1 



— = 7T 



nk^n'k' ,qa 



\vfaT 



[/ n v(*0-/»^)]«04'-fc-<r 

(18) 

Note that Eq. (fT8|) agrees with the Landau damping rate 
derived from the Fermi golden rule in the literature! 44 ! 45 
The above solution suggests that the Landau damping 
process can be understood as the resonant absorption of 
the SPP by electrons. The two <5-functions in Eq. ([TBI) in- 
dicate that for a monochromatic SPP wave with wavevec- 
tor q, the absorption occurs between pairs of states \nk) 
and \n'k'} satisfying the energy and momentum conser- 
vations 



nn 
^kk' 

k' -k 



0, 



(19) 
(20) 




(C) Low temperatures 

(strong) (weak) 

IPI 




6f=f k -f k 



High temperatures 




FIG. 2: (Color online) Schematic of the resonant pairs cor- 
responding to the case of (i) strong Landau damping and 
(ii) weak Landau damping in the electron spectrum (a) and 
in the electron distribution (b). The two wavevector re- 
gions of the resonant pairs are illustrated by thick green/blue 
curves corresponding to the strong/weak Landau damping 
regime. The resonant pairs are centralized around the res- 
onance corresponding to the SPP central wavevector Q s as 
indicated by the vertical black dotted curves. The horizon- 
tal dashed line marks the chemical potential of the electrons 
fi. (c) The electron polarization and electron population dif- 



ference as function of the center wavevector K 



k+k' 



with 



k' — k = Q s . (i)/(ii) corresponds to the strong/weak Landau 
damping regime at low temperatures and (iii) corresponds to 
both the strong and the weak Landau damping regimes at 
high temperatures. The vertical solid line marks the reso- 
nance corresponding to Q s . 
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Each pair of the states \nk) and \n'k') consist a resonant 
pair {nk,n'k') relevant for the SPP damping. For the 
multi-subband system, there usually exist several such 
resonant pairs, laying between different subbands n and 
n! and being well-separated from each other. For a non- 
monochromatic SPP wave packet with sufficiently narrow 
spectrum in q, the two states \nk) and \n'k') of each res- 
onant pair become two wavevector regions. Note that 
in the following discussion, we shall focus on the non- 
monochromatic SPP wave packet, and the resonant pair 
is referred to as the wavevector region unless otherwise 
specified. The resonant pairs are illustrated in Figs. HJa) 
and (b). By using the resonant pairs, one can rewrite 
Eq. (HI]) into 

r 1 = (21) 

i 

Tf 1 = ^E E \p^9f\ 2 [Mk')-fum 

q<? (nfc,n'fe')ei 

x «>*'-fc-<?= (22) 

with i being the index for the resonant pair corresponding 
to the SPP wave packet, whose spectrum is decided by 
the line-shape function p® s . (nk,n'k') 6 i means that 
the two states \nk) and \n k') belong to the i-th resonant 
pair. Thus one can see that the damping rate of the SPP 
wave packet is the sum over the absorption rates of all 
the relevant resonant pairs. 

Note that the electron population difference Sf = 
fn'a(k') — fna{k) of the corresponding resonant pair plays 
an important role on its contribution to the SPP damp- 
ing. For the degenerate electrons where a well-defined 
Fermi surface exists around the chemical potential, there 
exist two regimes of the SPP damping: (i) a strong Lan- 
dau damping regime where states \nk) and \nk') of a 
resonant pair lay in each side of the chemical potential in 
the electron spectrum, leading to a large population dif- 
ference Sf and hence a strong SPP damping; (ii) a weak 
Landau damping regime where the chemical potential lies 
outside all the resonant pairs, leading to a small Sf and a 
weak SPP damping. This is illustrated in Figs. 0(a) and 
(b). Note that at high enough temperatures, the Fermi 
surface can be smeared out and such difference vanishes. 

It should be emphasized that according to Eq. (TT5)) , the 
resonant pair can be visualized as the resonant peak in 
the polarizations between the two subbands n and n' . In 
Fig-HJc), we illustrate the polarizations P(nk,n'k') cor- 
responding to the resonant pairs shown in Figs. (Ha) and 
(b) as function of the center wavevector K = (k + k') jl 
with k' — k — Q s . One finds that the polarization exhibits 
a Lorentzian peak around the resonance corresponding to 
the central wavevector Q s , as indicated by Eq. (TT5l) . 

Note that such resonant peak can show different fea- 
tures in the strong and weak Landau damping regimes at 
low temperatures. In the strong Landau damping regime, 
the corresponding polarization P exhibits a strong reso- 
nant peak concentrated in the region with large Sf [(i) in 
Fig. He)], indicating a large SPP absorption by the elec- 



trons. In contrast, in the weak Landau damping regime, 
the corresponding resonant peak is weak and lies outside 
the large Sf region [(ii) in Fig. He)], indicating a small 
SPP absorption. In addition to the resonant peak, side 
peaks can also exist in the off-resonant regime due to the 
corresponding large Sf. At high temperatures where the 
Fermi surface is smeared out, the population difference is 
rather flat for both the strong and weak Landau damp- 
ing regimes and the corresponding polarization exhibits a 
strong peak around the resonance for both regimes [(iii) 
in Fig. He)]. 

D. Influence of the scattering 

The scattering influences the Landau damping by 
changing the resonant excitation of the polarizations. 
Specifically, (1) the scattering can introduce dissipative 
channels, inducing a decay of the polarization; (2) the 
scattering between polarizations with different preces- 
sion frequencies induces a frequency-mixing, leading to 
a modification of the polarization precession frequency. 
These two effects can be further clarified by assuming 
that for each resonant pair, the scattering term has the 
form 

Ink.n'k' =2ZT i [P a (nk-q,n'k'-q)-P <J {nk,n'k% (23) 

9 

where Ti stands for the phenomenological relaxation rate 
for the polarization of the i-th resonant pair. Note that 
(nfc, n'k') and (nfc— q, n'k' —q) belong to the i-th resonant 
pair. 

For each resonant pair, the polarization P in the scat- 
tering term I given above can be obtained by treating 
Ti perturbatively and solving Eq. (j!3[) order by order, 
yielding 

P a (nk,n'k';t) = &v- k - 9 gf ' p^B\{t)[f na {k) - /„,,(*')] 
x ^W'-It-HfJ)* - 1]/( W »»' - r? + if \), (24) 

where 

. .nn 

rf = ri£(i- **^), (25) 

q k—q,k'—q 

r? = ^ 5>$ # - (26) 

9 

The summation in Eq. (|25p is understood as a princi- 
pal value integral. Note that we have omitted the k, k'- 
dependence of f inside each resonant pair for simplic- 
ity. The detail of the derivation is given in Appendix [Cj 
By comparing Eq. ((24)) to Eq. (fT5|) , one can see that the 
detuning wPS in the resonant denominator is modified 
into wjgfc, —Tf, indicating that the polarization precession 
frequency is shifted by the scattering. The scattering 
also induces a finite imaginary part to the resonant 



6 



Low temperatures 
(strong) 




K 



(b) Low temperatures 

(weak) 



(j) broadening 



6f=f k -f k 




(C) 



High temperatures 
(strong)/(weak) 



(i) broadening 




(ii) 



shifting 



K 

(iii) 



K 



FIG. 3: (Color online) Schematic of the effect of the broad- 
ening and shifting on the electron polarization for (a) strong 
Landau damping regime at low temperatures; (b) weak Lan- 
dau damping regime at low temperatures and (c) both strong 
and weak Landau damping regimes at high temperatures. 
The thin solid vertical lines represent the resonance without 
broadening and shifting. The resonances with the broadening 
and shifting are represented by the thick solid and thin dotted 
vertical lines, respectively. For clarification, only the shifting 
towards the small K direction is illustrated. Note that the 
population differences Sf become flat in both regimes at high 
temperatures. 



denominator, representing the decay of the polarization 
due to the scattering. 

The above solution indicates that the scattering modi- 
fies the resonance between the polarization and the SPP 
by introducing both an energy shift and an energy broad- 
ening to the corresponding resonance pairs. On one hand, 



the energy shift modifies the energy conservation Eq. (fT9|) 
into 



J kk' 



r? = o. 



(27) 



Thus the corresponding resonant pairs are shifted by 
the scattering. On the other hand, the energy broad- 
ening loosens the energy conservation constraint given 
by Eq. (|19[) . Thus the corresponding resonant pairs are 
broadened. Note that the broadening and shifting are 
usually small and cannot induce overlaps between differ- 
ent resonant pairs. 

Accordingly, the broadening T\ and shifting f f also 
manifest themselves in the SPP damping rate. Following 
the same procedure of the derivation of Eq. (IT8|) . the 
SPP damping rate in the presence of the scattering can 
be written as 



(28) 



q& (nk,n' k')Gi 



W-rr) 2 + (f^ 5fe '- fc - 9 ' 



(29) 



By comparing the above equations to Eq. (|18[) . one ob- 
serves that the (^-function corresponding to the energy 
conservation Eq. (|19p is broadened into a Lorentzian with 
width and shift Tf. 

It should be emphasized that the broadening and shift- 
ing of the resonance pair can also be visualized as the 
broadening and shifting of the corresponding resonant 
peak in the polarizations as illustrated in Fig. [31 This of- 
fers a simply way to interpret the influence of the broad- 
ening and shifting on the SPP damping rate. The in- 
fluence of the shifting depends on the direction of the 
shift. From the figure, one can see that the shifting re- 
duces the resonant peak in the polarization if the peak is 
shifted towards the region with smaller Sf, thus the ab- 
sorption of the SPP by the corresponding resonant pairs 
is reduced [(iii) in Figs.^a) and (c)] and the correspond- 
ing SPP damping rate is suppressed. Otherwise, if the 
peak is shifted towards the region with larger Sf [(iii) 
in Fig. E)Jb)], the absorption is enhanced and the SPP 
damping rate is enhanced. 

The influence of the broadening can be different in the 
strong and weak Landau damping regimes at low tem- 
peratures as illustrated in Figs. |3fa) and (b). In the 
strong Landau damping regime, the broadening can re- 
duce the sharp resonance peak of the polarization [(ii) in 
Fig. [3fa)] and suppress the absorption of the SPP. Thus 
the corresponding SPP damping rate is suppressed in this 
regime. In contrast, in the weak Landau damping regime, 
the broadening of the resonance increases the absorption 
from the region with larger Sf [(ii) in Fig. [3f b)] , thus 
the absorption is enhanced and the corresponding SPP 
damping rate is enhanced. Note that at high tempera- 
tures, as the polarization exhibits a sharp peak around 
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the resonance in both the strong and weak Landau damp- 
ing regimes, the scattering tends to suppress the SPP 
damping rate in both regimes [(h) in Fig.[3jc)]. 



E. Broadening and shifting from a simplified model 

In order to gain a further understanding of the micro- 
scopic origin of the broadening and shifting, we discuss 
the contributions of the ei, ep, and ee scatterings to the 
broadening and shifting within a simplified model in this 
section. 

Within this model, we assume that the scattering only 
occurs between the polarizations inside each resonant 
pair. Under such assumption, all the scattering terms 
can be written in an unified form for the i-th resonant 
pair 



nk,n' k' 



^nk,n'k-M) p ^{nk,n'k')^ (30) 



where both (nk, n'k') and (nk — q, n'k' — q) belong to the 
i-th resonant pair. Note that the corresponding T a / b for 
each scattering mechanism can be obtained by comparing 
the scattering terms [Eqs. (|B1HB13)) ] to P° k n , k , given in 
the above equation. 

Note that Eq. ((30]) has a similar structure as Eq. ([23| . 
thus one can derive the corresponding broadening and 
shifting following the similar procedure, yielding 



V b = ( V b (a) — r a (n) 

i / ^ \ nk ,n' 'k' ,iVi J nk,n'k',i\H) 



J kk' 



, .Jilt 

W k-q,k'-q 



kk' 



k—q,k' —q 



)S(u 



k—q,k f —q 



(31) 



).(32) 



The broadening T b and shifting F" due to each scattering 
can be evaluated by substituting the corresponding scat- 
tering term into Eqs. (|30H32I) . Note that we have omitted 
the k , k' -dependence of inside each resonant pair 

to make the equation simple and physically transparent. 
For each resonant pair, T is evaluated by choosing 
(nk, n'k') corresponding to the SPP central wave vector 
Q s (i.e., k' — k = Q s ) since the line-shape function p® s 
is peaked at q = Q s . The shifting Tf can be written as 



-pa(ei) _ „ 
pa(cp) _ q 



^a(cc) 



2n x: E n «'K'°)' 

3=1,2 n 



(33) 
(34) 
(35) 



where 

n«,(nfc,g) 



with V q being the screened ee interaction. k\ 

-2(e«_ Qs - £ ™'+^)/g s -Q s /2and k 2 = 2(eJ-^'_ 0< 



n s )/Q s +Q s /2. f 



a(ci) -pa(op) 



and T ■ v ; correspond to the 



contributions from the ei, ep and ee scatterings, respec- 
tively. The corresponding broadenings have the forms 



f fe (o0 _ m * wr] . flQ ^0 \(„~,nn _ ~n'n'\ 

li - m TrnU |jfc/| ){v v Q ), 



(37) 



s6(ep) 



a Q 



X W£ofnAk -q) + N< Q f< a {k - q)]/\k - q\ 
Q 

X W^ofnAk -q) + N> }< a {k - q)]/\k - q\ 



q=q + 



q=ql 



+ {nk i — > n'k'}, 

= {2n^[U^nk,0)-KAnk,0)}/\k 



k\} 



nk 

+ {nk 4 — > n'k'}. 



(38) 



(39) 



where q { ± } satisfies [q { ± ] } 2 -2kq±±2m*n ho = and 5™"' 
represents the screened ei interaction, {nk < — > n'k'} 
stands for the same term as in the previous {} but with 
the interchange of indices nk < — > n'k' . 

From Eqs. (|33H39[) . one finds that different scattering 
has different contribution to the broadening and shift- 
ing. Only the ee scattering contributes to the shifting, 
while the contributions from the ei and ep scattering van- 
ish. Although all the scatterings can contribution to the 
broadening, their relative importance can be quite differ- 
ent. This is because for the nanowires considered here ; 
the ei, ep and ee matrix elements v™ 71 , M r Q q and V q nn 
are not sensitive to the subband index n. Thus, for the ei 
and ee scatterings, according to Eqs. (|37j) and (|39|) . the 
terms in the bracket can largely cancel each other, leading 
to small contributions left. In contrast, such cancellation 
is absent for the ep scattering, thus its contribution to 
the broadening is expected to be larger than the ones 
from the ei and ee scatterings. 

Equations (|28[) and (|2"9")l with the broadening and shift- 
ing given in Eqs. (|33H39p consist the analytic solution for 
the SPP damping rate. The analytic solution in this 
section offers a simple picture to understand the influ- 
ence of the scattering on the SPP damping: The SPP 
damping comes from the resonant absorption of the SPP 
by electrons, while the scattering can introduce a broad- 
ening and a shifting to the resonance and hence affects 
the damping process. At low temperatures, the broad- 
ening tends to suppress the SPP damping rate in the 
strong Landau damping regime. While in the weak Lan- 
dau damping regime, the broadening tends to enhance 
the SPP damping. At high temperatures, such difference 
vanishes and the broadening tends to suppress the damp- 
ing in both regimes. The shifting can suppress the SPP 
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damping if the resonance is shifted towards the region 
with smaller 5f, but boost it if the resonance is shifted 
toward the region with larger Sf. Different scatterings 
can have different contributions to the broadening and 
shifting. From the simplified model in this section, one 
finds that the shifting is determined by the ee scattering, 
and the broadening is mainly decided by the ep scattering 
for the typical nanowires considered here. 



III. NUMERICAL RESULTS 

In the numerical investigation, we choose nanowires 
to be free-standing InAs nanowires. The typical elec- 
tron density n is in the range of 10 17 ~ 10 18 cm -3 
and the wire radius R is around 25 ~ 75 nmjii The LO 
phonon energy J^lo = 29 meV and the electron effective 
mass m* = 0.023mo with mo representing the free elec- 
tron mass. The dielectric constants of the nanowires are 
ef 3 = 12.3 for high frequency and e\ = 15.5 for low fre- 
quency. The dielectric constant outside the nanowire is 
€2 = 1.0. Note that for such nanowires, there are 10-20 
electron subbands relevant to the SPP damping. We set 
the impurity line density = 0.5n e with n e = irR 2 fio 
being the electron line density. 

By numerical solving the kinetic equations Eqs. (fT2|) 
and ([TB")) . one obtains the temporal evolution of the 
SPP amplitude B s . The SPP damping rate r _1 can 
be extracted by fitting the real part of B s with a sin- 
gle exponential decay of the cosine oscillation: Re[B s ] = 
B° s exp(— tjr) cos(w s i), where the initial value of the SPP 
amplitude B® is chosen to be real. We set the wave packet 
length L = lOOi? (Ref.©. 



A. Landau damping: Size and temperature 
dependence 

Before we discuss the influence of the scattering, it 
is helpful to first obtain an understanding of the SPP 
damping without scattering. In Fig. Ufa), we show a 
typical behavior of the SPP damping rate as function of 
the SPP central wavevector Q s and wire radius R for the 
nanowire with electron density no = 1.5 x 10 17 cm" 3 . 
The temperature is chosen to be 100 K. One finds that 
the SPP damping rate oscillates with the radius R. Note 
that similar size-dependent oscillations have also been 
reported in metal nanoparticles and thin films i^ 3 - " 45 ' 76 " — 

Such oscillations are usually attributed to the quan- 
tized electron states in the nanostructures^ 3 -"— which 
can be understood in terms of the resonant pairs in the 
nanowires we studied here. To illustrate this, we con- 
centrate on the damping rate corresponding to a typical 
SPP central wavevector Q s = 1.66 x 10~ 2 nm -1 [skyblue 
curves in Fig. Eta)] and show the corresponding resonant 
pairs for radii R = 34, 38.5, 40 and 43 nm in Fig. 0|b). 
Each resonant pair can be represented by the resonance 
corresponding to the central wave vector Q s , since the 



line-shape function of the SPP wave packet is peaked at 
Q s . There are four resonant pairs in Fig. EJb) which lay 
between different subbands: pair (i) is between the sub- 
bands 1 and 4, pair (ii) is between the subbands 2 and 6, 
pair (iii) is between the subbands 3 and 9 and pair (iv) is 
between the subbands 4 and 8. In the figure, the four res- 
onant pairs (i-iv) are denoted by the skyblue dots, green 
squares, brown open circles and yellow triangles, respec- 
tively. The subbands corresponding to each resonant pair 
are also plotted with solid curves in the same color. Note 
that for R = 34 nm, only the resonant pair (i) is relevant 
for the damping. For R — 38.5 and 40 nm, both the 
resonant pairs (i) and (ii) are relevant. For R — 43 nm, 
all the four resonant pairs (i-iv) contribute to the SPP 
damping. 

From Fig. H^b), one can see that as the radius R in- 
creases, the resonant pairs move from left to right in the 
electron spectrum. When a resonant pair moves across 
the chemical potential marked by the horizontal black 
dashed lines, a crossover between the strong and weak 
Landau damping regimes occurs, which induces the size- 
dependent oscillations shown in Fig. Hta). Note that 
the peaks/valleys correspond to the strong/weak Lan- 
dau damping regimes. For example, the oscillation from 
R = 32 to 38.5 nm is due to the crossover induced by 
the resonant pair (i). While the crossover induced by the 
resonant pair (ii) induces the oscillation from R = 38.5 
to 43 nm. One finds from the figure that R = 34 and 
40 nm correspond to the strong Landau damping regime 
whereas R = 38.5 and 43 nm correspond to the weak 
Landau damping regime. 

One also observes from Fig.[4jb) that due to the many 
subbands in the nanowires, there usually exist multi- 
resonant pairs relevant for the SPP damping. For a given 
Q s , more and more resonant pairs become involved as 
the radius R increases. The number of relevant reso- 
nant pairs are labeled in the R-Q s plane in Fig. U{a). 
Note that as the number of the resonant pairs increases, 
the magnitude of the size-oscillations become less pro- 
nounced. This is mainly because the oscillations are usu- 
ally induced by the crossover due to one resonant pair as 
R varies. For the system with many resonant pairs, the 
contribution from one resonant pair becomes less signif- 
icant. Thus the size-dependent oscillations can be sup- 
pressed for nanowires with large R. 

It should be emphasized that the size-dependent os- 
cillations can also be suppressed by increasing temper- 
ature T. This is because the crossover is more pro- 
nounced for strongly degenerate electrons where a clear 
Fermi surface exists around the chemical potential. In 
high-temperature regime, the crossover is largely sup- 
pressed. To show this, we plot the damping rate r _1 
as function of the radius R and temperature T for Q s = 
1.66 x 10~ 2 nm -1 in Fig. HJc). One sees that as temper- 
ature increases, the damping rate r _1 corresponding to 
the strong Landau damping regime decreases, while r _1 
corresponding to the weak Landau damping regime in- 
creases. This leads to a suppression of the size-dependent 
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FIG. 4: (Color online) (a) The SPP damping rate as function of SPP central wave vector Q s and wire radius R without 
scattering for no = 1.5 x 10 17 cm -3 and T = 100 K. The skyblue curve represents the damping rate corresponding to the SPP 
central wave vector Q s = 1.66 x 10 -2 nm -1 . The number of the relevant resonant pairs are labelled in the R-Q B plane, (b) 
The resonant pairs for ii = 34, 38.5, 40 and 43 nm corresponding to Q s = 1.66 x 10 -2 nm -1 in the electron spectrum. For 
clarification, only the lowest 10 subbands are plotted, (c) The SPP damping rate without scattering as function of radius R 
and temperature T for Q s = 1.66 x 10~ 2 nm -1 . (d) Temperature dependence of the damping rate r -1 without scattering for 
R = 34, 38.5, 40 and 43 nm with Q s = 1.66 x 10 -2 nm -1 . Symbols correspond to the numerical results. Red curves represent 
the results from the analytic solution. The contributions of each resonant pair from the analytic solution are also plotted as 
curves with different colors. The skyblue double-dotted chain, green chain, brown dotted and yellow solid curves correspond 
to the contribution from the resonant pair (i-iv), respectively. 



oscillations in high-temperature regime. 

One can also obtain the above results from the ana- 
lytic solution of the kinetic equations without scattering 
[Eqs. (f2TT) and ([22])]. To show this, we compare the tem- 
perature dependence of r _1 from both the numerical (red 
squares) and analytic (red solid curves) solutions for the 
nanowires with ii = 34, 38.5, 40 and 43 nm in Fig. U{d). 
One finds good agreement between each other, indicat- 
ing that the analytic solution without scattering offers 
a good estimation to the numerical results. Note that 
according to the analytic solution, the temperature de- 
pendence of the SPP damping rate originates from the 
population difference of the resonant pairs. 

From the analytic solution, one can also identify con- 
tributions from different resonant pairs, which are plot- 
ted as curves with different colors and line shapes in 
Fig.Hfd). The skyblue double-dotted chain, green chain, 
brown dotted and yellow solid curves correspond to the 
contributions from the resonant pair (i-iv), respectively. 
It is clear that the relative importance of the resonant 



pairs can be quite different. For the strong Landau damp- 
ing regime, there usually exists one resonant pair whose 
contribution is much larger than the other pairs. For ex- 
ample, the damping rate t" 1 is mainly determined by the 
resonant pairs (i) and (ii) for ii = 34 nm and R = 40 nm, 
respectively. For the weak Landau damping regime, the 
contributions from different resonant pairs can be com- 
parable. For example, for R = 38.5 and ii = 43 nm, al- 
though the resonant pair (ii) has a large contribution to 
the damping rate r _1 , the other resonant pairs can also 
play important roles, especially at high temperatures. 



From the above results, one finds that the SPP damp- 
ing exhibits size-dependent oscillations and distinct tem- 
perature dependence without scattering, which can be 
explained by the analytic solution. 
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FIG. 5: (Color online) (a) SPP damping rate in the presence 
of all the scattering as function of radius R and temperature 
T for Q s = 1.66 x 10~ 2 am" 1 , (b) SPP damping rate r" 1 cal- 
culated from the numerical results for Q s — 1.66 x 10~ 2 nm -1 
with different scattering for R = 34, 38.5, 40 and 43 nm. Sym- 
bols with big blue dots, small green squares, small olive dots 
represent r _1 calculated with the ep, ee and ei scatterings, 
respectively. The brown triangles represent r _1 calculated in 
the presence of all the three scatterings and the red squares 
represent r^ 1 without scattering. 

B. Influence of scattering 

Now we discuss the influence of the scattering on the 
SPP damping. In Fig. [SJa), we plot the damping rate 
t -1 as function of the radius R and temperature T for 
Q s = 1.66 x 10~ 2 nm -1 in the presence of all the scatter- 
ing. Comparing to the case without scattering [Fig. [He)], 
one finds that the scattering has pronounced influence 
on the SPP damping: (1) The size-dependent oscilla- 
tions are effectively smeared out, and (2) the temper- 
ature dependence also becomes weaker compared to the 
case without scattering. 

To gain a better understanding of the influence of the 
scattering, in Fig.[5Jb), we compare the temperature de- 
pendence of the damping rate r _1 with and without scat- 
tering for nanowires with four typical radii R = 34, 38.5, 
40 and 43 nm. The brown triangles represent r _1 cal- 
culated in the presence of all the three scatterings, while 
t _1 without scattering are plotted with red squares for 
comparison. Note that R = 34 and 40 nm correspond to 
the strong Landau damping regime, whereas R = 38.5 
and 43 nm correspond to the weak one. 
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FIG. 6: (Color online) Comparison between the numerical 
and analytical results of the SPP damping rate for R = 34, 
38.5, 40 and 43 nm with the ep scattering. Blue dots repre- 
sent the numerical results, while the blue dashed curves show 
the results from the analytic solution. The contribution of 
each resonant pair from the analytic solution is also plotted. 
The skyblue double-dotted chain, green chain, brown dotted 
and yellow solid curves correspond to the contribution from 
the resonant pair (i-iv) , respectively. For comparison, the nu- 
merical and analytical results for the damping rate without 
scattering are also plotted with red squares and solid curves, 
respectively. 



In the presence of the scattering, it is seen that the 
damping rate t _1 is markedly suppressed in the strong 
Landau damping regime (R = 34 and 40 nm). In con- 
trast, for the weak Landau damping regime (R = 38.5 
and 43 nm), the scattering plays different roles in differ- 
ent temperature regimes: The damping rate is markedly 
enhanced in the low-temperature regime (T < 150 K), 
but is largely suppressed in the high-temperature regime 
(T > 150 K). A crossover exists at the intermediate tem- 
perature regime. It is also noted that the damping rate 
can be enhanced/suppressed by almost one order of mag- 
nitude by the scattering. 

To understand these influences, we first identify the 
dominant scattering mechanism. To do so, we calcu- 
late the damping rates r _1 with the ep, ee or ei scatter- 
ing only, and plot them with big blue dots, small green 
squares and small olive dots in Fig- E^b), respectively. It 
is clear to see from the figure that the damping rate r _1 
is dominated by the ep scattering^ 

We first concentrate on the ep scattering. From the 
analytic solution within the simplified model, we have 
attributed the effect of the ep scattering to the broaden- 
ing of the resonant pairs. To see if this picture gives a 
proper description of the influence of the ep scattering in 
general case, we calculate the temperature dependence 
of the damping rate r _1 by using the analytic solution 
Eqs. (|28p and (|2"5)) with the ep-scattering-induced broad- 
ening given in Eq. (|38|) . The calculated analytic results 
are compared to the numerical ones in Fig. 

In the figure, the blue dots represent the damping rate 
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t -1 from the numerical results, while the blue dashed 
curves represent t _1 from the analytic results. For com- 
parison, we also plot the numerical and analytic results 
without scattering with red squares and solid curves, 
respectively. One finds good agreement between each 
other, indicating that the broadening can give a proper 
description of the effect of the ep scattering. Note that 
at low temperatures, the broadening tends to suppress 
the SPP damping rate in the strong Landau damping 
regime (R = 34 and 40 nm). While in the weak Landau 
damping regime (R = 38.5 and 43 nm), the broadening 
tends to enhance the SPP damping rate. At high tem- 
peratures, such difference vanishes and the SPP damping 
rate is suppressed in both the strong and weak Landau 
damping regimes. This leads to a crossover between the 
suppression and enhancement for R = 38.5 and 43 nm 
corresponding to the weak Landau damping regime as 
shown in Fig. [TJ1 This also agrees with the conclusion 
from the analytic solution. 

It is worth noting that the scattering can also suppress 
the temperature dependence of the SPP damping rate. 
This is because the temperature dependence originates 
from the population difference of the resonant pairs. For 
the resonant pairs with the broadening, the correspond- 
ing population difference is less sensitive to the temper- 
ature, leading to the suppression of the temperature de- 
pendence of the corresponding SPP damping rate. 

We further point out that the scattering can also 
change the relative importance of different resonant pairs. 
To see this, we identify contributions of different resonant 
pairs from the analytic solution, as applied in Sec. IIII Al 
for the case without scattering. In Fig. |51 the skyblue 
double-dotted chain, green chain, brown dotted and yel- 
low solid curves correspond to the contributions from the 
resonant pairs (i-iv) , respectively. One finds that the res- 
onant pair (i) dominates the damping for R = 34 nm, 
whereas the resonant pair (ii) plays the most important 
role for R = 38.5 and 40 nm. For R — 43 nm, all the 
four resonant pairs have comparable contributions to the 
damping, with the largest contribution coming from the 
resonant pair (iii). Comparing to the case without scat- 
tering [Fig. SJd)], one also finds that the relative im- 
portance of different resonant pairs is modified by the 
broadening, especially in weak Landau damping regime. 

From the above discussion, one comes to the conclu- 
sion that the influence of the ep scattering on the SPP 
damping rate can be understood as the broadening of the 
resonant pairs. The analytic solution incorporating such 
broadening shows good agreement with the numerical re- 
sult. Note that in the above results, the contribution of 
the SO phonons is omitted since it is much smaller than 
that from the LO phonons. This is shown in detail in 
Appendix [DJ 

Now we briefly address the ee and ei scatterings. From 
the analytic solution within the simplified model, the ef- 
fect of the ei scattering is attributed to the broadening 
of the resonant pairs. While for the ee scattering, the 
main effect is due to the shifting. The SPP damping rate 
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FIG. 7: (Color online) Comparison between the numerical 
and analytical results of the SPP damping rate for the ee 
and ei scatterings for R = 34, 38.5, 40 and 43 nm. Small 
green squares and green dotted curves represent the numer- 
ical and analytical results for the ee scattering, respectively. 
Small olive dots and olive double-dotted chain represent the 
numerical and analytical results for the ei scattering, respec- 
tively. For comparison, the numerical and analytical results 
for the damping rate without scattering are also plotted with 
red squares and solid curves, respectively. 



with the ee/ei scattering can also be calculated from the 
analytic solution [Eqs. ([25)1 and (129)) with the broaden- 
ing and shifting given in Eqs. (|33)) and (l37l) for the ei 
scattering and Eqs. (|35|) and (|39| for the ee scattering, 
respectively]. The analytic results are compared to the 
numerical ones in Fig. [7J From the figure, one observes 
qualitatively good agreement between each other, indi- 
cating that the ee and ei scattering can also be under- 
stood as the broadening and/or shifting of the resonant 
pairs. 

It is pointed out that the effect of the broadening and 
shifting can be visualized from the behavior of the po- 
larization, which gives a more intuitive picture for the 
effect of the scattering. This is discussed in detail in 
Appendix [E] 

From the above results, one finds that the scatter- 
ing tends to smear out the size-dependent oscillations 
of the SPP damping rate. The temperature depen- 
dence can also be suppressed by the scattering. Note 
that this effect is quite general and can be seen for 
nanowires with different electron densities. To demon- 
strate this, we show the size and temperature depen- 
dence of the SPP damping rates for nanowires with elec- 
tron density no = 5.0 x 10 17 cm~ 3 without and with 
scattering in Figs. |8fa) and (b), respectively. The cen- 
tral wavevector of the SPP wave packet is chosen to be 
Q s = 1.3x 10~ 2 nm -1 . From the figure, it is seen that the 
size-dependent oscillations are smeared out by the scat- 
tering. The temperature dependence is also suppressed. 
These effects are similar to the ones for nanowires with 
no = 1.5 x 10 17 cm~ 3 investigated above, indicating that 
these effects are quite general for typical InAs nanowires. 
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FIG. 8: (Color online) Size and temperature dependence of 
the SPP damping rate for nanowires with electron density 
no = 5.0 x 10 17 cm -3 (a) without any scattering and (b) with 
all the scattering. The central wavevector of the SPP wave 
packet is chosen to be Q s = 1.3 x 10~ 2 am" 1 . 



onant pairs. Moreover, different scattering has different 
contribution to the broadening and shifting. The main 
effect of the ei/ep scattering is to cause a broadening, 
whereas the main effect of the ee scattering is to intro- 
duce a shifting. The effect of the broadening and shift- 
ing can be incorporated into an analytic solution, which 
shows good agreement with the numerical result. 

To demonstrate the effect of the scattering, we inves- 
tigate the damping of the axial symmetric SPP mode 
in InAs nanowires in the presence of the ei, ee and ep 
scatterings. Without any scattering, the SPP damping- 
exhibits size-dependent oscillations and a distinct tem- 
perature dependence. In the presence of the scattering, 
the size-dependent oscillations are markedly smeared out 
and the temperature dependence is also suppressed. The 
damping rate can be enhanced/suppressed by almost one 
order of magnitude. For InAs nanowires investigated 
here, the ep scattering is found to be dominant. These 
results are found to be general for typical InAs nanowires, 
which demonstrate the importance of the scattering on 
the SPP damping for semiconductor nanowires. It is 
further pointed out that our model can be applied to 
nanowires made of other semiconductors, offering a sys- 
tematic way to investigate the effect of electron scattering 
on the SPP damping in such systems. 
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IV. CONCLUSION AND DISCUSSION 

In conclusion, we present a microscopic kinetic theory 
to study the effect of the electron scattering on the Lan- 
dau damping of the SPP in semiconductor nanowires. 
Based on the semiclassical model Hamiltonian of the 
SPP-electron system and the nonequilibrium Green- 
function approach, we derive the kinetic equations of the 
SPP-electron system, with all the scattering explicitly 
included. Within this model, the SPP damping is un- 
derstood as the absorption of the SPP by the electron 
polarization of the resonant pairs. The population dif- 
ference of the resonant pairs Sf plays an important role 
on the SPP damping, leading to a strong and a weak 
Landau damping regimes for degenerate electrons. The 
scattering influences the SPP damping via the broaden- 
ing and shifting of the resonant pairs, which have dif- 
ferent effects on the strong and weak Landau damping 
regimes. At low temperatures, the broadening tends to 
suppress the SPP damping in the strong Landau damping 
regime. Whereas in the weak Landau damping regime, 
the broadening tends to enhance the SPP damping. At 
high temperatures, this difference tends to be vanished. 
The shifting can suppress the SPP damping if the reso- 
nance is shifted towards the region with smaller Sf, but 
boost it if the resonance is shifted toward the region with 
larger Sf. The broadening and shifting can be visualized 
from the corresponding electron polarization of the res- 
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Appendix A: Semiclassical model for SPP 

Our derivation of the quantized SPP Hamiltonian 
Hgpp and its coupling to electrons -ffspp-ci follows the 
procedure used in Rcf. 80, in which the SPPs are ob- 
tained from quantization of the classical plasmon field. 
The classical plasmon field is derived within the hydro- 
dynamical model, which treats plasmons as irrotational 
deformations of the conduction electrons From this 
model, the Hamiltonian of the plasmon can be written 

/m* 
dr—n{r)\Vr,{r)\ 2 

e 2 f e 2 f 
+ y / dr— / dr n(r)n(r)V ee (r,r) 

+ J drG[n(r)}, (Al) 

where the irrotational flow has been assumed, i.e., v(r) = 
Vr/(r) with v(r) being the velocity field. n(r) is the 
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electron density and to* represents the electron effective 
mass. 

For both r and r' inside the nanowire, the Coulomb 
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interaction V ee (r,r') has the form 
1 



V ee (r,r') = 



2 ,ei° 
r — r'\ ir £2 



e im(tp-tp') 



f 

JO 



dkcos[k(z - z')]C m (kR, —)I m {kp)I m (kp r 

C2 



(A2) 



where 



K m (kR)K' m {kR) 



(2 



I m {kR)K' m {kR) 



^I' m (kR)K m {kR) 
(A3) 

Note that the distortion of the Coulomb interac- 
tion V cc (r,r') by the dielectric mismatch between the 
nanowire and the environment has been taken into ac- 
count, which is not addressed in previous work. 

G[n(r)] represents the exchange, correlation and inter- 
nal energy of the electrons, which is approximated by the 
Thomas- Fermi functional 



G[n(r)} = 



10m' 



-(30 



2\2/3_5/3 



i) '"(r), 



(A4) 



with the exchange-correlation contributions neglected. 

From the above Hamiltonian, up to the first order of 
the perturbation, one can derive linearized hydrodynamic 
equations as 



d t n x {r) = V • [n (r)V7 7l (r)], 

e 

TO* 

57 ni(r) 



d t m( r ) 



dr'V ee (r,r') ni (r') 



3m* Mr)] 1 / 3 ' 



(A5) 



(A6) 



where n\ and 771 are perturbations around the equilib- 
rium. no(r) — hqQ(R — p) is the electron density in the 
equilibrium with hg being the average electron density. 

7 - 



5m* 



The normal modes from the above equations include 
both the surface plasmon and volume plasmon modes. 
The surface modes, which we focus on in this paper, can 
be represented by the ansata^ 3 . 



11qm{r) = Ry~]Qqm.e 



im.(p-\-iqz 



(qp) + AI m { K p)\, (A7) 



where I m (p) is the modified Bessel function of the first 
kind. A and k are parameters depending on to and q, 
with to being the angular quantum number of the SPP 
mode and q representing the SPP wave vector along the 
wire axis. Here the first term in the square bracket rep- 
resents the incompressible deformation of the electrons 
while the second term represents the correction due to 
the finite compressibility. 



1 



By substituting the ansatz Eq. (| A7|) into Eqs. (|A5 
and (|A6I) , one obtain the equation for the parameter k 



-?-(K 2 -q 2 )- (qR)I m+1 (qR)K m (qR) 

/<T' ( KRI m+l(KR) _ I m (KR) \ 
X i m \qRI m+1 (qR) I m (qR)J 
■ kR I m+ i{nR) 



-A(l + C' qm ) 
qR I m +i(qR) q 



/ KRI m+1 

V qR I m+ i 



(kR) I m (KR) K m+1 (qR) \ 



} 



(qR) I m+1 (qR) K m (qR) J\ 

0, (A8) 



with Lo p = 



- J 57 =2/3 



3m«"0 



.4: 



where 



pb- being the bulk plasma frequency and 
The parameter A can be expressed as 
I m +l(qR)(qR)(l + C qm ) 

J f^Bl/vmfl Xf (qR)I m+ i(qR)I m (K,R.) \ ' 

(A9) 



1 



1 + 



e? Im+l(qR) K m (qR) 
C2 I m (qR) K m + 1 (qR) 



(A10) 



Once k is obtained, the corresponding eigen-frequency 
£lqm can be calculated from the equation 



1 + /3 2 



"<jm 



J 'qm 



(KR)I m+1 (KR)I m (qR) - (qR)I m ( K R)I m+1 {qR)] } 

= 0. (All) 



Note that without the dielectric mismatch (e^° = £2), 
Eqs. (|A8IIAll|) agree with the results for the surface mode 
from Ref.il. 

The Hamiltonian of the SPP can be obtained by sub- 
stituting Eq. (|A7[) into Eq. (|A1[) and keeping only terms 
up to the linear order, which can be written as 



H 



2TTR d m*n 



where the canonical coordinate Q' qm has the form 

Qm 



' qm 



AI m (nR)n + I m (qR)q 



AI m (nR) + I m {qR) Q q 



(A13) 



Following the canonical quantization procedure, the col- 



lective coordinate Q' qm can be transformed into 



■> qm 



27rR 3 m*noQ q7 



qm. 



— qm) 1 



(A14) 
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where the annihilation(creation) operator b qm (b qm ) sat- 
isfies the boson commutation relation [b qm ,b qlm ,] — 
5q, q >5 mm >. Substituting Eq. (|A14|) into Eq. (|A12|) . the 
Hamiltonian of the SPP is quantized as 



Eq. (J5)), the corresponding matrix element v q n can be 
calculated as 



e 2 S~ ~, 



1 Jrn(\tp)Jrh,(\%'p) 



" 2^e?Rj Q ^Ji h+1 (\t +1 )J fn , +1 (\f; +1 ) 

Hspp = ^!Vi(2 fe t m 5 9m + 1) . (A i5) x[l (qRp)K (qR) + C' (qR)I (qRp)K (qR)]. (A21) 



By choosing m = 0, one obtains the Hamiltonian of the 
axial symmetric SPP mode under investigation 



(A16) 



without taking into account the zero-point energy which 
is irrelevant for our study. The angular quantum number 
m is also omitted without confusion. 

Following the similar procedure, the induced potential 
of the SPP mode can be quantized as 



V m (p^,z) 
where 

K in (p,^) = 



dqV q in ( Pl ^z)(b q + bl q ), (A17) 



"o 



y 2Trm*uj p ^T q 
I (qp)C' q ' + AI (k P ) 



^/[Ah{nR)nR + I 1 (qR)qR][AI ( K R) + I Q (qR)} 

(A18) 



with 



C" - - 



A 



{kR)AI\ (kR) + (qR)h(qR) 
K Q (qR)+C' Q (qR)K Q (qR) 
(KR)I 1 (KR)K (qR) + {qR)I (nR)K 1 (qR) 



C' (qR)K (qR)((KR)Ii(KR) - (qR)I (KR) 



h(qR) Y 
Io(qR)J- 
(A19) 



The SPP-electron coupling Hamiltonian ifspp-ci in 
Eq. © can be obtained by combining the single electron 
states in Eq. ([5]) and the induced potential of the SPP 
given above. The corresponding matrix element reads 



9k-k> = J pdpd^dz^{p,L Pl z)Vl n _ k ,{p^,z)i} n < w {p^,z) 1 

(A20) 

with ipnk given in Eq. (|5|). 

Now let us discuss the matrix elements for the inter- 
action Hamiltonians for electrons. For the ei interaction, 
we have assumed that the impurities are distributed on 
the surface of the nanowire with an axially symmetric 
distribution. 74 Given the electron wavefunction ip n k in 



Note that the dielectric mismatch effect has been taken 
into account in the above expression. For the ep inter- 
action, we here consider the polar interaction between 
the electrons and the LO phonons. The corresponding 
matrix element reads 

V27re 2 f2 L o(l/ef°-l/e;) 2 



M nn = 

Qg J (\m\ T ( \m' 



x / dppJ A (Xtp)J^(Xf;p)e l ^ R 6 Arh ,, (A22) 



where q — (Q, q) represents the LO phonon wave vector 
and f^LO is the LO phonon energy. e p is the unit vector 
along the radial direction of the wire and is the dielec- 
tric constant of the nanowire in the static limit. Note 
that here we model the LO phonons by 3D modes, which 
is adequate for nanowires with large diameteri 68 ' 69 For 
the ee Coulomb interaction, the corresponding matrix el- 
ement reads 



2e 2 



i ,i 
dpipi / dp 2 p2 



IoiqRp^KoiqRpy) 



+ l)I (qRpi)Io(qRp 2 )Co(qR, — ) 



</m + l(A™) Jm' + l(A™ ) 



(A23) 



where p> = max(pi,p2) and p < — min(pi, /c^)- Here 
we have omitted the interband term of the Coulomb in- 
teraction which can be small for nanostructures^ Note 
that the dielectric mismatch effect has been taken into 
account in the ee Coulomb interaction. 



Appendix B: Derivation of kinetic equations 

The time evolution of B q (t) — (b q (t)) can be derived 
from the Heisenberg equation of the annihilation oepra- 
tor b q , which reads 

d t B q (t) = -iQ q B q (t)+ p2°- k ,9k-k>G<{n'k',nk-,tt). 



nn' .kk' .a 



(Bl) 



By using the definitions of B s and P given below Eq. (|14[) , 
one has 



d t B s (t) = -iJ2(n Qs+q -n s )p^B Qs+q e 



kk' 
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where the first term in the right hand side of the equation 
describes the effect of the SPP dispersion on the dynam- 
ics of the SPP wave packet and the second term describes 
the dissipative effect due to the coupling to electrons. For 
the SPP wave packet with sufficiently large length L, the 
line-shape function p® s exhibits a sharp symmetric peak 
around the central wave vector Q s . In this case, one can 
perform the Taylor expansion of the dispersive relation 
Q,Q a+q around Q s up to the linear order of q to elimi- 
nate the first term. In doing so, we omit the effect of 
the SPP dispersion on the wave packet dynamics and 
consider only the dissipative effect. 

Within the nonequilibrium Green-function approach, 
the quantum kinetic equation of electrons can be derived 
following the Kadanoff-Baym method , 66 i 67 yielding 

- id t - (e# - ejj)l G| (nk, n'k'; tt) 

= Y.( B -i + B D [5r'G|(nfc,nfc' - q;tt) 

nq 

-g™Gf(nk + q,n'k';tt)} + I% n , k ,(t), (B3) 

where the first term in the right hand side of the equation 
is the coherent driving term of electrons by the SPP and 
the second term describes the scattering 



G${kk;tt)Y>{kk' ;tt)) 



X 



}, (B4) 



where [>< — y<] stands for the same term as in the pre- 
vious [ ] but with the interchange of >< — X. It should 
be emphasized that due to the driving of the SPP, the 
dressed Green function G§ has off-diagonal components 
with respect to the momentum k and subband index n. 

By treating the driving term as a perturbation, the 
kinetic equation can be linearized by keeping only terms 
up to the linear order, 



id t -(el, - el) G^(nk,n'k';tt) 

= (B k _ k , + Bl_ k )[g^: k Gl(nk,nk;tt) 

- 9^- k Gl(n'k', n'k'; tt)] + J** (t). (B5) 



The first term in the right hand side of the equation is the 
linearized driving term. Note that within the lineariza- 
tion, the dressed Green function G& in the driving term is 
replaced by the free electron Green function Gq ct , which 
is diagonal with respect to the momentum k and sub- 
band index n. By using the definition P a (nk,n'k';t) = 
~iG<{nk,n'k';tt)e inst and f na (k) = -iG^(nk,nk;tt), 
one obtains Eq. ([TU)) from the above equation. 

Now we turn to the scattering term. We first 
discuss the ei scattering. Following the standard 
approximation,— the corresponding scattering term can 



be expressed as 



nk.n' k' 



S ci , a (nk, n'k'; >, <) - Si a {n'k', nk; <, >) 



(B6) 



S ci , a {nk,n'k' ;>,<) = J dt^^mv^v^' 
~°° nn'n kq 

x G>(nk - q,n'k - q;tt)G<(nk,n'k';tt), (B7) 



where v q n 



y nn j e nn ^ j g screenec j electron- 

impurity interaction. The screening e nn (q) is evaluated 
following Ref.[83 in the static limit, which can be written 

as 



y q n V ^fUk+ q )-fUk) 



la 



'k+q 



i0+ 



(B8) 



Note that S e - h<7 contains the product of two dressed Green 
functions. By keeping only terms up to the linear order, 
S cl a can be linearized as 



~nn 
Q 



S ei , a (nk,n'k'\>,<) = [ dt^^mv^v, 
J -°° fin - kg 
x GQ a (nk — q,nk~ q;tt)G< (hk, n' k! ;it) 



J ~°° aa> r 



nn' k q 

x G>(nk - q, n'k' - q; tt)G^ a {n 'k' ', n'k'; tt). (B9) 

Note that within the linearization, S C i ia is separated into 
two terms. Each term contains only one dressed Green 
function Ga , while the other one is replaced by the free 
Green function G^ a . By applying the rotating wave 
approximation with respect to the SPP frequency fi s , 
one can remove the non-frequency-matched terms in the 
above equation, yielding 



,(nk, n'k'; >, <)e 



~in s t 



nn'q 



n s ) 



«f<™ fUk-q)P.{n'k,n'k') 

~ v™v* n> f< a {k')P„{nk - q, n'k' - q)] (BIO) 



in the Markovian limit. Substituting Eq. (|B10|) into 
Eqs. (|B6|) and (|B7|) . and using the definition Inhn'k' = 



nk.n' k 



, h ,e 



, one gets the final expression of the impurity 



scattering term 



nk.n' k' 



£ {*& 



n s ) 



nn'q 



" P °( nk ~ 9' n ' k ' - ?) - vf ^T Pa{n'k, n'k') 
+ 5(st:_ q -e n k -n s ) 



P a {nk - q, n'k' - q) - V™ ' v n q n ' P a {nk, nk')] } 



(Bll) 
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The ep scattering term can be derived following the sim- 
ilar procedure, which reads 



= * E E - 4 + n - + 

nn' qQ 



X 



M£M&' (N> f<Jk>) + N< f> a (k') 
x P a {nk + q, nk' + q) 

- M$pl™ (N> f> a (k + q) + N< f < a (k + q) 

x P a (n'k,n'k') 

+ 5{el +q - si +n s - n LO ) 

x 'M%M$? (iV< /< CT (fc') + N> Q f> a (k')) 
x P a (nk + q, n'k' + q) 

- M%M$ (N< f> a (k + q) + N> f<„(k + q)) 

x P a (n'k,n'k') 

(N> f< a (k) + N< Q f> a {k)) 
x P a (nk + q, nk' + q) 

'Q'i - u Qv 

x P a (nk, nk' 

+ 5(st: +q -si-n s -n LO ) 

'm%M%? (iV< /„< (k) + N> Q f>„(k) 
x P a {nk + q, n'k 1 + q) 
tf'Mg (N< f> a (k' + q) + N> f< a (k' + q) 

x P a {nk,nk') }, (B12) 



MA'/ MA*' (N> f> a {k' + q) + N< f< a (k' + q) 



-M, 



where N^ >] = N LO + \ - (+)± with N LO = 
l/[exp(^4^) — 1] representing the thermal LO phonon 
distribution. The ee scattering term can also be derived 



in the similar way, 

qnk 

x 'f<, a (k')Ul,(nk,q)P a (nk-q,n'k' - q) 
-f> a (k-q)n^(nk,q)P a (nk,n'k' 

+ 5{el_ q -ei+n s +el-el_ q ) 
f> l<T {k')W£{hk, ~q)P a (nk - q, n'k' - q) 
- f< a (k - q)nZ(nk : -q)P a (nk,n'k') 
+ 5{el~ei_ q + n s +el-el +q ) 
f< a {k)Xl n nl (nk,q)P a {nk - q, n'k' - q) 
~ fn'Ak' - q)WS(nk,q)P r7 (nk,n'k') 

f> a ( W(nfc, -q)P a {nk - q, n'k' - q) 

~ fn>.(k' ~ q)^(nk, -q)P a {nk, n'k')] }, 

(B13) 

where IL^(nk,q) = £ CT V™V™' f^{k + q)fa(k) with 
ynn _ ynn j ^nn ^ being the screened ee interaction. 

Appendix C: Scattering-induced frequency 
modification and decay of the polarization 

The polarization P for the i-th resonant pair in the 
presence of the scattering term I° k given in Eq. (|23|) 
can be solved by treating 1^ as perturbation and solve 
Eq. (|13p order by order. By assuming 

oo 

P a (nk,n'k') = J2 P a j) (nk,n'k';t), (CI) 

j=o 

one obtains 

d t P^ {nk, n'k') = ujffi'pW ( nk > n ' k '; *) 
+ i9k-k'P^ k 'Bt [fna(k) - f > (C2) 
for 0-th order and 

d t P ( J +1) {nk,n'k') = iul$P^ +1 \nk,n'k';t) 
+ J2 ^[P^Hnk - q, n'k' - q) - P^(nk, n'k')], (C3) 

for (j + l)-th order, with j > 0. 

In the long time limit t — > oo, the solution of the above 
equations reads 

P^(nk,n'k') = -gT-' k <ptk>Bl 

x [f n *(k) - kWW' + *0+), (C4) 
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for 0-th order, and 

P^ +1 )(nk,n'k') = P - ] ^ n ' k ') 



._ ^ r p£* ) (n*-«,» / * / -?) 

*1 j / 7T7 

Y Py\nk,n'k') 



1], (C5) 



for (j + l)-th order. 

Now let's evaluate J^ g [- 
order, 



(nk—q,n'k'—q) 
1?!p {nk,n'k') 



11. For 0-th 



n 



Pi°\nk- q,n'k' - q) 
P!?\nk,rik') 

gr'k>Pklk>Bnfn«(k -q)- /„v(fc' - q)} 



l.: I.:' 



k—q,k—q' 



,nn 
J kk' 



k—q,k—q f 



i0+ 



(C6) 



where we have assumed that f n <j(k) varies slowly with 
respect to k and can be cancelled out. The 1-st order 
term can be evaluated as 



P^\nk- q,n'k' - q) 
P^\nk,n'k') 



- 1 



- (nfc - g, n'tf - g)]{^[Pi 0) (nk - q, n'k' - q) 

Q 

nn' . 

-pV>(nk,n'k')]}- 1 lj 



k—q,k—q' 



iO 



(C7) 



where we have assumed that ^2 q [Pa°\nk — q, n'k' — q) — 

Pa (nk, n'k')] varies slowly with respect to k(k') and can 
be cancelled out. 

Following the similar procedure, for j-th order, one has 



E 



pj j) (nk - q, n'k' - q) 



p( j \nk,n'k') 



Ei 



k—q,k—q' 



-1] 
1], 



(C8) 



with j > 1. 

By combining Eqs. ([OT> . (IC4l) . (IC5l) . (IC6l) and (IC8l) . 
one obtains Eqs. (|24ti26[) . 
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FIG. 9: Temperature dependence of SPP damping rate in 
the presence of the ep scattering for R — 34 and 38.5 nm. 
Blue dashed curves represent the results with only the LO 
phonon scattering, while the brown solid curves represent the 
results with both the LO and SO phonon scatterings. The 
damping rates without any scattering are plotted with red 
double-dotted chain curves for comparison. Other parameters 
are all the same as Fig. [6] 



Appendix D: Effect of the SO phonons 

The SO phonons in nanowires§2 can also influence 
the damping of the SPP. However, for the typical InAs 
nanowires we considered here, the SO phonons are of 
marginal importance since their contribution is much 
smaller than that from the LO phonons. This will be 
shown in this Appendix. 

The induced potential due to the SO phonon is given 

b}i 88-90 

y so = ^(T^e^e^a^ + h.c), (Dl) 
p 

with 



^so 



/27re 2 P D(Cl 



so) 



i P (qr), 



D(u) 
£i( w ) 



y q I p (qR)P p (qR) p 

62 M 

£2(^)^761(0;) - ei (u)d-e 2 (uy 



O 2 

S2 Ti 



or 



(D2) 
(D3) 
(D4) 



where Op° is the annihilation operator for the SO 
phonons and f2so is the SO phonon energy. fiLi(^Ti) 
represents the LO(TO) phonon energy inside (i = 1) 
and/or outside (i = 2) the nanowire. The coordinates 
are chosen according to Fig. [TJa). 

As the SO phonon energy is very close to the LO 
phonon energy for typical InAs nanowires^ we assume 
ilso = ^lo in the calculation. By using the single elec- 
tron states in Eq. ([5]), the corresponding electron-SO- 
phonon interaction Hamiltonian can be written as 



Pso — 



pq nn'ka 



[a s ° + (al 



p-q 



)^] C nka C n'k- 



q<j 



,(D5) 
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where the interaction matrix element reads 

pq ~ J^+i(Af)J^ +1 (Af/) 6rhSn ' +p 

x f dppJ^Xt 9)^(^9) , Ip{qRf>) ■ (D6) 
Jo ^I p (qR)I p (qR) 

By using the analytic solution given in Sees. Ill Dl 
and III El we calculate the temperature dependence of the 
SPP damping rates for R = 34 and 38.5 nm in the pres- 
ence of the ep scattering with and without the contribu- 
tion of the SO phonons, which is plotted in Fig. [9j One 
can see that the damping rates are dominated by the LO 
phonons and the SO phonons have very little effect. 



Appendix E: Visualization of the broadening and 
shifting 

We have seen that the SPP damping rate can be un- 
derstood as the broadening and shifting of the resonant 
pairs. In the analytic solution, the broadening and shift- 
ing can be seen clearly from both the resonant denomi- 
nator in the polarization [Eq. (|24l) ] and the Lorentzian in 
the damping rate [Eq. ([29]) ] . One may wonder whether 
the broadening and shifting can also be observed in a 
more intuitive way from the numerical results. In fact, 
inspired by Eq. (|2"4"1) . the broadening and shifting of the 
resonant pairs can be visualized from the normalized po- 
larization P(nk,n'k')/(Sf ■ B s ), from which the structure 
of the resonant denominator can be identified. 

Before we show the numerical results of the normal- 
ized polarizations, we first explain what we expect from 
the polarizations. According to the analytic solution 
Eq. (|15p . the magnitude of the normalized polarization 
\P(nk,n'k')/(5f ■ B s )\ without scattering can exhibit a 
sharp Lorentzian peak around the resonance correspond- 
ing to the SPP central wavevector Q s as indicated by the 
resonant denominator. Several side peaks can also exist 
as the SPP wave packet is nonmonochromatic [e.g., the 
line-shape function pj s ^ 6(q — Q s )}- With scattering, 
these peaks can be shifted and broadened, indicating the 
shifting and broadening of the corresponding resonant 
pairs. Small side peaks may also be smeared out by the 
broadening. Note that with the scattering, the polariza- 
tions from the numerical results can have more complex 
behaviors. The peaks in the polarization can also be dis- 
torted by the scattering, which has been omitted in the 
analytic solution. 

We first concentrate on the typical case with R = 
34 nm and T = 150 K corresponding to the strong Lan- 
dau damping regime. The magnitudes of the normalized 
electron polarization \P(nk, n'k')/(Sf-B s ) \ from the com- 
putation are plotted in Fig. [TUJa). In the figure, polar- 
izations with different scatterings are plotted by curves 
with different colors. The contour plots of the polariza- 
tions are also shown in the K-q plane, which are useful 




FIG. 10: (Color online) 3D plot and the corresponding con- 
tour plot of the normalized (a) and unnormalized (b) elec- 
tron polarization corresponding to the resonant pair (i) for 
R = 34 nm, T = 150 K. The corresponding population differ- 
ence is also plotted in the figure. The olive, green and blue 
curves represent the polarization with the ei, ee and ep scat- 
tering only, respectively. The red curves represent the polar- 
ization without scattering, while the brown curves represent 
the polarization in the presence of all the three scatterings. 
In the contour plots, the resonance corresponding to the SPP 
central wavevector Q s is shown with skyblue dot in the con- 
tour plot. The population difference is plotted with orange 
curves. Different curves are offset along g-axis for clarify. To 
provide a clear visualization, in (a) the normalized polariza- 
tions with the ep and all the scatterings are enlarged by a 
factor of 2. In (b), the unnormalized polarization with the 
ee/ep/all scattering is enlarged by a factor of 4/30/30. 



for identifying the shape and position of the polariza- 
tions. The population differences Sf = f n '{k') — f n (k) for 
the resonant pairs are also plotted with orange curves for 
comparison, with the corresponding contour map plot- 
ted in the K-q plane. The polarizations with different 
scatterings and the population difference have been off- 
set along the g-axis for clarification. The corresponding 
unnormalized electron polarizations \P(nk, n'k')/B s \ are 
also plotted in the similar way in Fig. UOf b) for compar- 
ison. Note that all the polarizations are taken at time 
t = 2.86 ps. Polarizations taken at other time show sim- 
ilar behaviors. 

In Fig. [TUKa) , the polarization is plotted for the reso- 
nant pair (i) which is the only relevant resonant pair for 
the damping. The resonant pair is centralized around the 
resonance corresponding to the SPP central wavevector 
Q s , which is represented by the skyblue dots in Fig. rTOT a). 
One finds that without scattering, the corresponding po- 
larization (red curve) exhibits a sharp main peak around 
the resonance. Several side peaks exist around the main 
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FIG. 11: (Color online) 3D plot and the corresponding con- 
tour plot of the normalized electron polarization correspond- 
ing to the resonant pair (ii) for (a) R — 38.5 nm, T = 100 K, 
(b) R = 38.5 nm, T = 150 K and (c) R = 38.5 nm, T = 200 K. 
The corresponding population difference is also plotted in 
the figure. The resonance corresponding to the SPP central 
wavevector Q s is shown with green square in the contour plot. 

peak. These features agree with the previous discussion. 
Note that the main peak lies in the regime with large 5f, 
indicating that the SPP absorption by the polarization 
is strong. Also note that due to the strong resonance, 
the profile of the unnormalizcd polarization is mainly de- 
cided by the resonance and shares a similar structure as 
the normalized one. 

In the presence of the ei scattering, the polarization 
is slightly modified. Some side peaks are smeared out. 
However, the broadening is rather small and the main 
peak is almost unchanged. This can be better seen by 
comparing the polarization with the ei scattering (olive 
curve) to the one without scattering (red curve). Ac- 
cording to the previous discussion, these features indicate 



that the ei scattering introduces a small broadening. 

In the presence of the ee scattering, a shifting of the 
peaks can be seen clearly by comparing the contour plot 
of the corresponding polarization (green curve) to the 
one without scattering (red curve). Note that small 
side peaks can also be identified with the ce scattering. 
These features indicate that the ee scattering introduces 
a large shifting to the corresponding resonant pairs. The 
broadening due to the ee scattering is small since the 
small side peaks are not smeared out. It is also noted 
that the ee scattering can distort the polarization. The 
peak becomes fragmented and the profile becomes non- 
Lorentzian. However, as the influence to the SPP damp- 
ing comes from the summation of the polarizations as 
indicated by Eq. (fT2")l . the effect of these distortion on 
the SPP damping is marginal. 

In the presence of the ep scattering, the peak in the 
polarization (blue curve) is markedly broadened and all 
the side peaks are smeared out, indicating that the ef- 
fect of the ep scattering introduces a large broadening. 
Note that the maximum of the broadened peak is also 
shifted compared to the one without scattering. This 
shows that the ep scattering can also have contribution 
to the shifting of the resonant pairs. However, due to the 
large broadening, the effect of such shifting on the SPP 
damping is marginal. We also point out that due to the 
large broadening, the profile of the unnormalized polar- 
ization is mainly determined by the population difference 
rather than the resonance, which can be seen from the 
corresponding unnormalized polarization (blue curve) in 
Fig. Hlb). 

In the presence of all the three scatterings, one finds 
that the polarization (brown curve) has almost the same 
profile as the one with the ep scattering only. The broad- 
ening and shifting due to the ei and ee scatterings are 
negligible due to the large broadening introduced by the 
ep scattering. This indicates that the ep scattering plays 
the dominant role. 

For the weak Landau damping regime, as the scat- 
tering has different influence on the SPP damping rate 
in different temperature regimes. We plot the normal- 
ized polarizations in Fig. [IT] for three typical cases: (a) 
R = 38.5 nm, T = 100 K, (b) R = 38.5 nm, T = 150 K 
and (c) R = 38.5 nm, T — 200 K, corresponding to the 
low, intermediate, and high temperature regimes, respec- 
tively. The polarizations are plotted for the resonant pair 
(ii) which dominates the SPP damping in these cases. 
Similar broadening and shifting can also be identified in 
the figure. Note that in these cases, the shifting due to 
the ee scattering is rather small and difficult to be iden- 
tified in the figures. 
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